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Abstract 



Optical images of the nearby star Fomalhaut show a ring of dust orbiting the central star. This dust is in many respects expected to be 
similar to the zodiacal dust in the solar system. The ring displays a clear brightness asymmetry, attributed to asymmetric scattering of 
the central starlight by the circumstellar dust grains. Recent measurements show that the bright side of the Fomalhaut ring is oriented 
away from us. This implies that the grains in this system scatter most of the light in the backward direction, in sharp contrast to the 
forward-scattering nature of the grains in the solar system. In this letter, we show that grains considerably larger than those dominating 
the solar system zodiacal dust cloud provide a natural explanation for the apparent backward scattering behavior. In fact, we see the 
phases of the dust grains in the same way as we can observe the phases of the Moon and other large solar system bodies. We outline 
how the theory of the scattering behavior of planetesimals can be used to explain the Fomalhaut dust properties. This indicates that 
the Fomalhaut dust ring is dominated by very large grains. The material orbiting Fomalhaut, which is at the transition between dust 
and planetesimals, can, with respect to their optical behavior, best be described as micro-asteroids. 
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1. Introduction 

Optical images of some nearby stars show the presence of 
orbiting solid material through the starlight it scatters (e.g 
iKalas et al., 2005). This has been interpreted as evidence of an 
evolving planetary system. Collisions between large bodies such 
as asteroids are believed to be responsible for the production 
of these small solid particles. This is similar to our own so- 
lar system, in which interplanetary dust grains are found that 
are thought to originate from asteroids and comets. The size 
and chemical composition of these grains determine the way in 
which they reflect starlight. Laboratory measurements and light 
scattering theory show that small solid particles scatter most in- 
cident light in the forward direction, in agreement with astro- 
nomical observations. However, recent observations of the dust 
orbiting the nearby star Fomalhaut have shown for the first time 
that the grains in this system scatter most of the light in the 
backward direction (Le Bouauin 'et all l2009l) . This is in sharp 
contrast with the current understanding of light scattering and 
challenges the analogy between interplanetary dust and the dust 
observed in Fomalhaut. In this paper, we show that an excellent 
match to the observed backscattering behavior of the material or- 
biting Fomalhaut can be obtained by using the theory of regolith 
covered surfaces. 



2. Apparent backward scattering 

The circumstellar ring of Fomalhaut as observed by the Hubble 
space telescope cl early shows a dark side and a bright side 
(IKalas et al.L l2005 i Regarding the orientation of the ring, the 
general forward scattering behavior of dust grains has led to 



the conclusion that the bright side of the ring should be the 
one tilted towards the observer. Using the Very Large Telescope 
Interferometer (VLTI), the rotation axis of the central star and 
subsequently the i nclina tion axis of the system was inferred 
dLe Bouquin et al.L 12009 ). This inclination is opposite to pre- 
vious assumption in that the direction of scattering is reversed 
from forward-dominated to backward-dominated. This contrasts 
with the current consensus in the field of light scattering. 

A relatively simple explanation for this can be found by con- 
sidering bodies much larger than interplanetary dust grains. As 
an example let us take the Moon, which for all practical pur- 
poses can be considered a gigantic dust grain. However, it is ev- 
ident that the Moon is bright when viewed in a backward scat- 
tering situation (full Moon) and it is dark in the forward scat- 
tering case (new Moon). The question then arises: is the Moon 
truly backscattering? The answer is "no". The apparent surprise 
that the Moon is forward scattering is caused by the sharp and 
in most cases invisible diffraction peak it causes. For the Moon 
this diffraction spike, caused by interference effects due to the 
shadow of the Moon, is so narrowly forward peaked that it can 
only be seen in the forward-most fraction of a degre^l i.e. when 
the Sun is right behind the Moon. For all practical purposes, the 
Moon appears to be backscattering, and diffraction can be omit- 
ted. 

For smaller grains, the diffraction spike often dominates the 
overall scattering behavior over a broader range of scattering an- 
gles. It appears that the disk around Fomalhaut contains dust 
grains that behave backward-reflecting like asteroidal bodies. 
This must be because the grains are very large, making the 
diffraction spike so narrowly forward-peaked that it is outside 
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Figure 1. Angular scattering functions for the regolith particles, 
showing the relative intensity scattered at a scattering angle 9. 
The pure reflectance function (solid, black curve) is given, and 
diffraction effects are added for particles with a narrow size dis- 
tribution around 3 /im grains (dashed curve). The exact scatter- 
ing function for this narrow size distribution of spheres as ob- 
tained from the Mie theory is also shown (solid, grey curve). 

the observable range of angles. For the Fomalhaut system, which 
has an inclination of 25°, the range o f obser vable scattering an- 
gles is 25° < < 155°. lKalasetal.1 (120051) already noted that 
the apparent single scattering albedo of the dust grains is in the 
range of 5 - 10%, which is low for any kind of dust grains. Such 
values are more typical for the geometric albedo of asteroids. 
The geometric albedo of asteroids does not include the contri- 
bution from diffraction. For large grains, it is easily shown that 
when diffraction is incl uded, the single sca ttering albedo is al- 
ways larger than 50% dvan de Hulstni957L Chapter 12). Thus, 
the true single scattering albedo of the grains, averaged over all 
angles and including diffraction, must be much higher than the 
observed value, indicating that most of the scattered light is not 
detected. 

3. Modeling the scattering function of asteroidal 
bodies 

Below we treat the dust grains in the Fomalhaut system as 
micro-asteroids, meaning that we use the theory of reflectance 
by regolith-covered asteroidal surfaces to compute their scatter- 
ing behavior For the reflectance of asteroidal surf aces, we use a n 
analytical model for the bidirectional reflectance (Hapkel. [l98lh . 
For simplicity we use the theoryj3mitting effects of macroscopic 
roughness as discussed in iHapkd ( Il984[) . This basically means 
that we set the macroscopic roughness parameter as discussed 
in that paper to zero. In addition to this, we ignore the oppo- 
sition eff'ect. This mainly plays a role near backward-scattering 
and thus largely falls outside the range of scattering angles that 
we have access to in the Fomalhaut system. To integrate over the 
surface of the body, we use Monte Carlo ray tracing. 

3.1. The regolith particles 

As input in these computations, we need to know the albedo and 
angular scattering function for the regolith particles that cover 
the surface of the grain. Note that since the regolith particles are 
closely packed, we need to k now these par ameters omitting the 
contribution from diffraction (lHapkelll981h . The contribution of 
diffraction can be separated when we compute the angular scat- 
tering function and single scattering albedo using geometrical 
optics. More exact computational methods like the Mie theory 
do not allow the separation of the contribution from diffraction. 
We thus assume the regolith particles are large enough so that ge- 



ometrical optics can be applied, and in addition we only consider 
surface reflections, i.e. transmitted rays are ignored. We realize 
that using geometrical optics is a rough approximation for the 
grain sizes we are considering. Nevertheless, since we are inter- 
ested only in the overall shape of the angular scattering function 
mostly at intermediate scattering angles, the approximation is 
not so bad for grain sizes a few times the wavelength of radia- 
tion, especially for particles with a con s iderab le imaginary part 
of the refractive index (Wi elaard et al.L |1997|) . In this approxi- 
mation the scattering properties of the particles are independent 
of their size and shape as long as they have a convex surface. To 
compute the reflections we need to specify the refractive index 
of the regolith particles which we take to be m - 1.6 + 0.1/, 
thought to be typical for cosmic materials. However, we find 
that the final result is hardly influenced by modest variations 
in the refractive index (with the real part ~ 1.5-1.8 and the 
imaginary part < 0.3). The angular scattering function for the 
regolith particles which results from this is dominantly forward- 
scattering, while their single scattering albedo is 1 1 % (omitting 
diffraction). In Fig. [T] we show three different curves: the angu- 
lar scattering function omitting diffraction (solid, black curve), 
the angular scattering function with added contribution from 
diffraction by a narrow size distribution of spheres with diam- 
eters around 3 fj.m (dashed curve), and the full angular scattering 
function computed using the Mie theory for this narrow size dis- 
tribution (solid, grey curve). All curves are plotted as a function 
of the scattering angle, i.e. 0° is forward-scattering while 180° is 
backward-scattering. Note that the curve obtained using reflec- 
tion plus diffraction is quite similar to that obtained using the full 
Mie theory. It can be argued that 3 fim is rather large for the re- 
golith particles covering the grains we discuss below. However, 
it does allow the grains to be covered by or composed of these re- 
golith particles. The parameters derived above are used as input 
for the optical properties of the regolith particles in the Hapke 
theory. 

3.2. The grains as a whole 

The resulting single scattering properties of the regolith particles 
obtained in the previous paragraph are inserted into the Hapke 
theory. Note that we use the theory in its most simple form, i.e. 
we ignore effects of macroscopic roughness and the opposition 
effect. In that case, the only remaining parameters are the single 
scattering albedo and the angular scattering function of the re- 
golith particles obtained above. As we will show below, we can 
find an almost perfect fit without considering additional parame- 
ters on macroscopic roughness and opposition effect. This means 
that using only the measurements available so far, we cannot 
constrain them and decided to use the most simple form that can 
still reproduce the observations. Much more accurate observa- 
tions of the angular scattering function of the Fomalhaut grains 
are needed to compare all Hapke parameters to those derived for, 
for example, the Moon or solar system asteroids. 

In Fig. |2] we show the resulting angular scattering function 
of the model Fomalhaut dust grains. The curves show the rel- 
ative intensity scattered at an angle 6. The grey area indicates 
the observable range of scattering angles for the Fomalhaut sys- 
tem. Note that the grains we simulate using the Hapke theory ap- 
pear to be predominantly backward-scattering (we see the cres- 
cents of the grains), while the small regolith particles covering 
the surface of these larger grains are predominantly forward- 
scattering, as discussed above. The angular scattering function 
resulting from only using the Hapke reflectance theory, i.e. with- 
out including diffraction, is shown by the long dashed line. To 
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Figure 2. Angular scattering functions for large dust grains, showing the relative intensity scattered at a scattering angle 6. The pure 
reflectance function (long-dashed curve) is given, and diffraction effects are added for 10, 30 and 100 / im grains (dotted, dashed 
and solid curves). The empirical scattering function of the Fomalhaut disk grains from Kala s et alj ( I2005I) is also shown (thick 
grey curve, albedo chosen to match the computations). The grey area indicates the observable range of scattering angles for the 
Fomalhaut system. 



simulate the finite size of the dust grains, we added the con- 
tribution from diffraction using a diameter of the grains of 10, 
30 and 100//m in dotted, dashed and solid lines, respectively. 
For the diffraction, we took simple Fraunhofer diffraction by 
a spherical aperture. To get rid of the resonance structures as- 
sociated with Fraunhofer diffraction, averaging was performed 
over a narrow size distribution (a flat distribution from 0.8 to 
1.2 times the average size). Grains of 100/im diameter clearly 
fit the best to the empirical angular scattering function (shown 
by the thick gray line) in agreement with Dent et al. (2000). For 
the observed scattering function of the grains orbiting Fomalhaut 
we used th e Heny ey-Greenstein parameterization obtained by 
iKalas et al.l (l2005l) . where we simply switched the forward and 
backward scattering directions, i.e. switching the asymmetry pa- 
rameter g from -1-0.2 to -0.2. The angular scattering function ob- 
tained in this way was multiplicatively scaled to match the com- 
putations. This scaling factor gives the single scattering albedo 
of the grains. We found that when we integrated the angular scat- 
tering function resulting from the Hapke theory over the angles 
accessible in the Fomalhaut system (the area indicated in gray), 
the albedo of the grains averaged over the observable range of 
angles is 5%. Thi s is consistent with the low albedo found by 
IKalas etal.l(l2005h . 

If we consider the assumptions on the regolith particles and 
the Hapke theory it is clear that the model we present is quite 
rough. However, a firm outcome of the model is that for the range 
of scattering angles we observe the scattering is dominated by re- 
flection rather than by diffraction. This implies, as shown above, 
that the grains are larger than lOOjum in diameter. This outcome 
does not depend on the assumptions in the model and is therefore 
a solid conclusion. 

4. Conclusion 

We conclude that the scattering surface in the disk around 
FormaUiaut is dominated by grains of at least 100 yum in size. Our 
findi ngs agree with the fact t hat the infrared spectrum is feature- 
less dStapelfeldt et al.Ll2004i) . indicating that grains smaller than 



a few micron are heavily depleted in the system. Also.'Pe nt et al.l 
(2000) find that large grains dominate the thermal emission. It is 
remarkable that such large grains are observed directly in opti- 
cal light in an astronomical object. The reason for this is that 
while the dust mass in many systems is dominated by large 
grains (see e.g. iWilner et all l2005l : [Testi et al.L l2003h . they ai-e 
normally over-shone by a component of small grains. This com- 
ponent, though less massive, dominates the optical cross section, 
and therefore appearance, of the disk. The fact that very large 
grains are so visible on a global scale in the Formalhaut disk 
means that small grains have been cleaned out from this disk ex- 
tremely efficiently, and that all there is left are these very large 
grains, directly probed by reflected stellar light. 

We conclude that we are observing the transition part of pa- 
rameter space, going from dust scattering to planetesimal reflec- 
tion. In fact, one might conclude that, almost 400 years after 
Galileo observed the crescent of Venus, we are seeing the cres- 
cents of the large dust grains in the disk around Fomalhaut. 
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